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ABSTRACT 

We study the form factor f+ (H = D, B) involved in H-n transi- 
tions. We use data on D° — > tt~1 + v to find at low recoil. Then 
we use perturbative QCD methods to calculate at larger recoil the 
contributions to f? that break heavy quark symmetry (HQS). Us- 
ing HQS relations we obtain an estimate of f+ which includes first 
order corrections in l/(2mg). Comparison with recent data on 
B° -► n'l+u gives V uh = 0.0025 ± 0.0005 ei ± 0.0004 t/l 



1. Introduction. 



Heavy quark effective theory (HQET) provides a model independent framework to 
describe the hadronic interactions of heavy mesons |]J. In the infinite quark mass (ttiq) limit, 
HQET makes manifest extra symmetries of the QCD Lagrangian: a flavor symmetry relating 
matrix elements of mesons with different heavy quark content, and a spin symmetry that 
relates amplitudes involving mesons with the same heavy quark content but with different 
heavy quark spin. Corrections breaking the symmetry are parametrized as an expansion in 
1/(2i71q), and will be small if the momenta exchanged between the heavy quark and the light 
degrees of freedom inside the meson are smaller than tuq. 

In particular, the semileptonic decays H — ► 7r / v (with H = B,D) have been studied 
at first order in l/(2mg). For vanishing lepton mass the decay rate depends only on the form 
factor f^, 

(7r(p) I qYQ I H(v)) = (m H v» + jf) + f H (m ff / - p") , (1) 

where v is the heavy meson four velocity and p the pion momentum. In these decays the 
kinematical variable v-p will vary from m n (for a final pion at rest in the rest frame of the 
initial heavy meson) to m#/2 + m^/(2m#) (the highest recoil). At first order in l/(2m,Q), 
f+ can be expressed JJ] 

f H = V^H , J_ B) (2) 

J+ v-p + A H yM 2m Q n y ! 

where Ah = tuh* — ttih and cq m is a renormalization factor. A and B above are functions of 
v-p that do not depend on the mass of the heavy quark: they will describe D—n and B-tt 
transitions at the same recoil point v-p. In the infinite mq limit only A contributes to f? 
and one has f+/ f+ = \JmBlmn H . 

At low v-p the momenta exchanged between the heavy quark and the light degrees 
of freedom during the transition are of order m# — niq = A m Aqqd- I n consequence, the 
corrections breaking heavy quark symmetry (HQS) will be small, with B/A = 0(A). At 
large v-p, however, it will be necessary to exchange large momenta to keep the expectator 
quark inside the outgoing pion. These hard proceses will be able to distinguish the mass of 
the heavy quark and will introduce symmetric and non-symmetric contributions of the same 
order: B/A = 0(v-p) = O(toq). Since large v-p is the dominant kinematical region in heavy 
meson decays, it will not be possible to give there an accurate relation of with fP based 
on HQS only. Our objective in this letter is to determine the value of the contributions that 
break HQS and grow (relative to the symmetric contributions) with v-p. This will allow a 
relation between the two form factors valid at order B/A = 0(A) also in the region of large 
v-p. 



The calculation of f+ has been attempted using different methods. Soft pion relations 
based on partial conservation of the axial vector current imply single pole behaviour at low 
v-p |. This means 

fH _ g _ Cm H /2 



where q 2 = (p B ~ Px) 2 is given in terms of v-p as q 2 = m 2 B + m 2 — 2m B v-p and C is the value of 
the form factor at q 2 = 0. In the soft pion limit || single pole dominance holds at first order 
in l/(2mg): dominant terms and first order corrections follow separately this dependence on 
v-p at low recoil. QCD sum rules || support the single pole behaviour in the whole region of 
momentum transfer in B and D to tt decays, and light-cone sum rules || suggest that order 
l/(2mg) corrections are small, around a 10% in D decays. Analogous results, consistent with 
single pole dominance, are obtained || in the combined limits of chiral symmetry and large 
N c . All these arguments would imply that the functions A and B in Eq. (fj) are just constants 
at v-p < m B /2. On the other hand, the calculation of f+ using perturbative QCD methods 
(see 0-0) gives qualitatively different results. In B decays the values obtained seem to be 
smaller than expected 0, and the HQS-breaking contributions are dominant in most of the 
phase space. 

We pretend here to obtain from . This will require the evaluation of non- 
symmetric corrections (contributions to B/A in Eq. (§)) that grow with v-p and become 
of order 1 at v-p = O{mo). At large v-p these corrections involve the exchange of large 
momenta, which suggests that they can be analyzed using perturbative QCD methods. Then 
our strategy will be the following. First we will use data on D° — > n~l + v to normalize the 
total form factor at low recoil. Then a perturbative QCD calculation (we will take an error 
of order AA from soft contributions) will give us B at large v-p. From that we will obtain f? 
for all the values v-p accesible in B decays. This result will be used to evaluate the partial 
decay rate r(£>° — ► it~l + v), 

dr(ff° -> n-l+u) _ G 2 F m B 2 2 2 s B 2 

d(vp) ~ 12tt3 lKbl [[VP} mJ ' (4) 

and comparison with recent data will give us a value for the CKM mixing V u b. 

2. Form factor at low recoil. 

To find the form factor we use BR(D° ^l+u) = (3.8 ± 1.2) x 1(T 3 0. The 
differential decay rate of this mode can be read from Eq. (2) just by replacing m B — > and 
Vub — > Vus- Taking a mean life td = (0.415 ± 0.004) x 10~ 12 s and single pole dominance from 
v-p = m n to v-p = m D /2 + m^/(2mo) we obtain 



cL) 13.8^/BR(D0 -> tt-I+u) (0.85 ± 0.12) 



t D = v i = ^ ' (ft) 

J+ vp + 0.142 v-p + 0.142 ' v ; 



As discussed before, the deviations from single pole behaviour at the recoils in D decays are 
small, and can be neglected respect to the quoted experimental uncertainty. 



3. Perturbative QCD calculation. 

In the Brodsky-Lepage formalism [ IT | the heavy-to-light transition amplitude is ex- 
pressed 

(n(p) | q^Q | H(v)) = J dx dy Tr^(y,p) T^(vp,x,y) ^ H (x,m H v)} , (6) 

where (i = H, it) is an effective quark- ant iquark distribution amplitude (x and y are 
the fraction of longitudinal momentum carried by the quarks), q and Q = c,b c are light and 
heavy valence quarks, T£ is a hard scattering amplitude given at lowest order as the sum 
of collinear skeleton graphs in Fig. 1, and Tr stands for a trace over spin, flavour and color 
indices. Contributions from Fock states with extra qq pairs will be suppressed by powers of 
l/v-p. The tensor structure of is 

%(x,p) = ^~ <Mx) K + 0(A) + p] 75 , (7) 

where Ic is the identity in color space and the distribution amplitudes satisfy 

dx <f>i{x) = ^=/i (8) 

(with this normalization the pion decay constant is f n = 132 MeV). The sum of the two 
diagrams in Fig. 1 gives 

(n(p) | urQ I H(v)) =-Jdidy g 2 s Tr(T a T a ) ^-f n ^^-f H Tr[ 

— S 
(eA- $)i^{p'q + m Q )^ a {m H + m H ^757" + 

u g u Q 



D 9 D q 1 

(9) 

where now (p i are normalized to one, all the 0(A) uncertainties in the mass term of the meson 
wave functions are contained in the 0(1) factor e, r Ti{T a T a ) = 4, and we have written the 
fraction of momentum carried by the light quark in the heavy meson as x = D g and 

D q ^Q are the denominators in the gluon and quark propagators (see below). 

At the energies involved in B decays the amplitude above has the tendency to depend 
strongly on the treatment of infrared regions. A selfconsistent calculation will require taking 



into account double-log (Sudakov) corrections fl2]| , which exponentiate and suppress contri- 



butions where an internal quark or gluon is near the mass shell. Two infrared regions need 
to be considered. 

(i) In Eq. ([]) the dependence of and Tff on the transverse momenta k± of the partons 
relative to the parent meson have been neglected. However, this is not a good approximation 
near the end points of x and y, where the gluon proopagator becomes singular. Using a 
modified factorization formalism it has been shown [TI| that as the momentum transfer grows 
the configurations with a large (in a 1/Aq C d scale) quark-antiquark transverse distance b are 
suppressed by Sudakov corrections and do not contribute to the transition amplitude. Small b 
corresponds essentially to large k±, its Fourier transform variable. When x, y — > these values 
of k± are the typical momenta carried by the exchanged gluon in the scattering amplitude 
and provide an adecuate renormalization and factorization scale. The denominators in the 
gluon and ligth quark propagators, 

D g = -2A(x^.p + ^|) = -2AA g 



D q = -2A(xvp+^) = -2AA q , 



(10) 



will be regulated at the end points by a generic transverse momentum (acting as an infrared 
cut off) k 2 ± = (1.2-1.6 Aqcd) 2 PI- 



(ii) The internal heavy quark in Fig. 1(a) carries a momentum pq = p\ — y p2 + ki- 
Then, its propagator will be singular for a certain distribution of parton momenta in the 
outgoing pion: 

_ A 2 + k 2 

D Q = p 2 Q -m 2 Q = -2m H (yvp-A + — i) (11) 

It is clear that for y m A/ (v ■ p) the heavy quark will be near the mass shell even if the 
transverse quark-antiquark distance is forced by Sudakov corrections (and also by the explicit 
k\ distribution ||14|| ) to be small. However, The heavy-to-light vertex will have double- log 



corrections, generated by the diagram in Fig. 2(b), that will suppress this infrared region as 
well. The external light quark in Fig. 1(a) has an offshellness A q w Aqcd, then this correction 
will suppress the configurations where the heavy (internal) quark is also near the mass shell. 
For = v-p and S = \q±\/niQ we obtain a Sudakov factor 

S ~ exp [ - (In* 2 -pf. - ^■p{yv.p-A) ) } ^ 

6 71 " h-QCD A-QCD 

with fi = y v-p — A. S provides in the diagram in Fig. 1(a) an effective cut-off of the region 
in y with (y v-p — A) < Aqcd- 



We can now proceed with the perturbative calculation. In the expressions above we 
can separate HQS symmetric and nonsymmetric contributions. We write the decay constant 
and the heavy quark propagator as 

1 1^1 A 2 + k\, 



(l-TT-^T^) 



D Q 2m H A Q 2m Q A c 

(13) 

with Aq = y v-p — A. Substituting in Eq. (||) we obtain the HQS-breaking contribution to 



AA g 

5 A 2 + fc 2 _ _ 

{ ^~ ( c / X — ~( 4 ^ - 4A - % e A + 2e m*)] + 

— [c f (-x + 2e) + 4xv-p~ 2^- - 4x e A] } , 



(14) 



where /i 2 = 2A x y v -p + k\ 



To evaluate B we take A QCD = 0.2 GeV, A = 0.4-0.5 GeV, e =0-1, c f = -1.45 and 
(/hV^h) (0) =0.34-0.42 GeV 3 / 2 (deduced from / B = 160-200 MeV and {f B y/rnS)/{fDy/rn^) = 
1.4 ± 0.1 [|l^]). For the 7r meson we shall use the asymptotic distribution amplitude <j)^(y) = 
6y(l — y) [fL6~l| . Due to the inclusion of k\ , the result does not depend on the details of the 
distribution amplitude of longitudinal momentum in the heavy meson. For example, when 
the fraction x of momentum carried by the light quark (in units) varies between 0.5 and 2, 
at v-p = me/2 we have that B goes from 0.40 to 0.15. We shall asume a constant distribution 
between these two values of x and zero otherwise. At v-p = m^/2 the 50% of the result comes 
from contributions with a s < 0.5. 

In Fig. 3 we plot our perturbative result for B (dots) for v-p > ms/2. We observe 
that at these intermediate and large v-p B is constant, which means that contributions to f? 
that break HQS follow a single pole behaviour l/(v-p). This is also the behaviour of B given 
by the soft contributions at low |. We shall use the perturbative value to normalize B 
at v-p = rriB/2 and will add a constant soft contribution B = AA as an error. To estimate A 
we will use the experimental form factor deduced from D° — > -k~1 + i> data. We extrapolate to 
lower recoils the value B = 0.28 GeV 3 / 2 obtained perturbatively and substract its contribution 



to the total form factor in Eq. (6). Then we obtain A = 0.52 GeV 1 / 2 . This value of A will be 
constant as far us single pole dominance is a good approximation (v-p ~ ms/2 according to 
QCD sum rules @, §). 



The errors in the perturbative calculation of B will come from the variation of (A, k 



a) (27%), from the error in (/ff\/^ff) (10%), and from higher twist and higher order in a s 
corrections (25%). There is also a 8% error from order (I/itlq) 2 corrections generated in the 
diagrams in Fig. 1 (the corrections appearing at higher twist will be smaller, since they will be 
suppressed by powers of a s ). Soft contributions introduce an uncertainty AB = AA = 0.23 
GeV 3 / 2 (82%). Adding the errors in quadrature we obtain our low-recoil estimate of B: 



B = (0.28 ± 0.25 th ) GeV 3/2 . 



(15) 



At large v-p (3> me/2) the soft contributions would drop and B would converge to its 
perturbative value B = (0.28 ± 0.09' h ) GeV 3 / 2 . 

Since A has been deduced from the difference between the experimental form factor 
and the calculated HQS-breaking correction, its uncertainty will come from the error in B (a 
18%) and from the experimental error in fP (a 18%): 



A = (0.52 ± 0.09 ex ± 0.09 th ) GeV 1/2 , 



(16) 



4. f? and determination of V u t,. 



We can now estimate fP: 



fi 



vp + A B 



with p = 25/3. At v-p < m B /2 



a s {m c ) 
a s (m b ) 



2/0 



(A + ^-B) 



f 



;i.39±0.24 e:c ±0.25 
vp + 0.046 



th\ 



(17) 



Finally, using recent data on B° — > n l + v from the CLEO Collaboration |T7j we 
can estimate the mixing angle V u {,. For BR(5° — > 7t~1 + v) = (1.8 ± 0.5) x 10 -4 and r^o = 
(1.56 ± 0.06) x 10" 12 s we obtain 



V„ h = 0.0025 ± 0.0005 ex ± 0.0004 



iii 



(19) 



where the experimental error comes from uncertainties in the D and B to 7T I v branching 
ratios. 



5. Conclusions. 



We have analyzed the form factor f? (H = D, B) describing H-7r transitions. In the 
semileptonic decays of a heavy meson the momentum transfer is of order v-p ~ mg/2, so a 
priori the corrections breaking the HQS are of order 1. We have computed perturbatively these 
corrections using the Brodsky-Lepage formalism and including in our calculation Sudakov 
effects, which are necessary to obtain a selfconsistent result at these values of v-p. 

We normalize the form factor A at low recoil using data on D° — ► 7r~/ + z/, and we 
normalize the contributions B breaking HQS using a perturbative calculation. At intermediate 
and large v-p these contributions follow a single pole dependence, which is also the behaviour 
followed by soft contributions at lower recoils |J. We introduce in B an error of order AA 
from soft contributions that does not increase with v-p. 

Our results for ff and the partial decay rate T(B° — > tt~1 + i>) are slightly higher than 
other estimates. For example, 



/^( g 2 = 0) = 0.26 [18]; 0.33 0; 0.27 |20 



and 



r(s c 



TT-I+V) 



x 10 13 s" 1 



0.51 [18]; 0.74 0; 1.0 



(20) 



(21) 



versus our values 0.52 and 1.8, respectively. We should stress, however, that our values are 
normalized by the experimental data on D° — > 7r~Z + z/. In this sense, it may be significant 
that typical calculations of the analogous form factor and partial decay rate in D decays fall 
short respect the observed values: 



f"(q 2 = 0) = 0.5 [18]; 0.69 0.67 



(22) 



and 



T(D° -> n~l + u) x 10 10 s" 1 = 0.39 [18]; 0.68 ||; 0.80 @ (23) 

versus 0.91 and 0.92 from the experiment (where single pole dominance has been assumed to 
deduce the value of the form factor). 

Since we have actually computed the difference between and /+ given at first order 
in the mass of the heavy quark by HQS, the total values of the form factors that we deduce 
are strongly correlated. It is possible to redo the analysis in terms of the branching ratios 
and express the mixing V uh as 



V, 



ub 



BR(B° 



7T 



-l+u) 



(97 ± 16) JBR(D° -> ir-\+v) - (0.57 ± 0.50) 



(24) 



For the central experimental values we obtain V u b = 0.0025, with a 18% error from uncertain- 
ties in our calculation. 
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Figure 1: Heavy-to-7r hard scattering amplitudes. 




Figure 2: Tree-level heavy-to-light vertex and one- loop diagram giving the double- log cor- 
rection. See PI for a definition of light-cone variables and Feynman rules. 
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Figure 3: Form factor B. Dots correspond to the perturbative values obtained with our 
calculation. 



